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We assess when electron trapping nonlinearities are expected to be important in Langmuir waves. The basic 
criterion is that the effective lifetime, td, of resonant electrons in the trapping region of velocity space must 
exceed the period of trapped motion for deeply-trapped electrons, t b = (n e /Sn) 1 ^ 2 2Tr/u} pe . A unitless figure 
of merit, the "bounce number" Nb = td/rs, encapsulates this condition and allows an effective threshold 
amplitude for which Nb = 1 to be defined. The lifetime is found for convective loss (transverse and longi- 
tudinal) out of a spatially finite Langmuir wave. Simulations of driven waves with a finite transverse profile, 
using the 2D-2V Vlasov code LOKI, show trapping nonlinearity increases continuously with N B for side loss, 
and are significant for N B ss 1. The lifetime due to Coulomb collisions (both electron-electron and electron- 
ion) is also found, with pitch-angle scattering and parallel drag and diffusion being treated in a unified way. 
A simple way to combine convective and collisional detrapping is given. Application to typical underdense 
plasma conditions in inertial confinement fusion targets is presented. The results show convective side loss is 
usually the more potent detrapping process in a single F = 8 laser speckle. 
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I. INTRODUCTION 

The nonlinear behavior of Langmuir waves (LWs) is 
a much-studied problem in basic plasma physics from 
the 1950s to the present. In this paper, we focus 
on nonlinearity due to electron trapping in the LW 
potential well. This intrinsically kinetic effect has 
motivated theoretical work such as nonlinear equilib- 
rium (Bernstein-Greene-Kruskal (BGK) modes 1 ), Lan- 
dau damping reduction 2 -, nonlinear frequency shifter— , 
and the sideband instability^. Important applications 
of trapping occur in LWs driven by coherent (e.g., laser) 
light, including the laser plasma accelerator— and stim- 
ulated Raman scattering •: SKS;'* ! 1 . The latter allows 
the prospect of laser pulse compression to ultra-high am- 
plitudes (the backward Raman amplifier) 12 . In addi- 
tion, SRS is an important risk to ICF— both due 
to loss of laser energy and the production of energetic 
(or "hot") electrons that can pre- heat the fuel. Ignition 
experiments at the National Ignition Facility (NIF)3^ 
have shown substantial Stimulated Raman backscatter 
(SRBS) from the inner cones of laser beamsi&. The 
current study is prompted primarily by SRS-driven 
LWs. Much recent work has focused on nonlinear ki- 
netic aspects of SRS, including "inflation" due to Lan- 
dau damping reductioniZr— , saturation by sideband 
instability 22 , and LW self- focusing in multi-D particle- in- 
cell simulations^—, Vlasov simulations^, and theory 2 ^. 
One goal is to find reduced descriptions, such as enve- 
lope equations, that approximately incorporate kinetic 
effects2&2&. 
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Our aim is to provide theoretical estimates for when 
electron trapping nonlinearities are important in LW dy- 
namics. These allow for self-consistency checks - or inval- 
idations - of linear calculations of LW amplitudes. This 
work is therefore not intended to augment the litera- 
ture on nonlinear LW dynamics, although we do present 
Vlasov simulations to quantify the onset of trapping in 
the presence of convective transverse loss. We consider a 
single, monochromatic wave with electron number den- 
sity fluctuation 5ncos(kx — ut), and unitless amplitude 
SN = Sn/n e where n e is the background electron density. 
Let "resonant electrons" denote those moving sufficiently 
near the phase velocity v p = uj/k that they can poten- 
tially execute trapped (as opposed to passing) orbits. 
Deeply-trapped electrons have an (angular) frequency 
lub = uJpeSN 1 /' 2 (u}p e — n e e 2 /eom e defines the plasma fre- 
quency in SI units) , known as the bounce frequency, cor- 
responding to a bounce period tb = 2tt/ujb oc SN^ 1 / 2 . 
For a process that detraps electrons in characteristic time 
td, the unitless "bounce number" Nb = tdl T B measures 
how many bounce orbits a resonant electron completes 
before being detrapped. 

Our estimates stem from the assumption that nonlin- 
ear trapping effects are significant when Nb is roughly 
unity. Trapping nonlinearities develop continuously with 
wave amplitude, and are not an instability with a 
hard threshold. Vlasov simulations presented below of 
driven LWs with a finite transverse profile demonstrate 
this. In addition, comparison with transit-time damp- 
ing calculations 31 show the reduction in Landau damp- 
ing varies continuously with Nb and obtains a 2x reduc- 
tion for Nb ~ 1. Bounce number estimates are some- 
what qualitative and demonstrate basic parameter scal- 
ings. The quantitative role of trapping depends on the 
specific application. 
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We consider here two detrapping processes: convective 
loss and Coulomb collisions. For a LW of finite spatial 
extent, electrons enter and leave the wave from the sur- 
rounding plasma (assumed here to be in thermal equilib- 
rium, i.e. Maxwellian) . Trapping will only be effective if 
these electrons complete a bounce orbit before transiting 
the wave. We find the detrapping time for both longitu- 
dinal end loss, which can be important in finite-domain 
ID kinetic simulations, and for transverse side loss in 2D 
and 3D. To quantify the effect of trapping in a LW with 
finite transverse extent, we perform 2D-2V simulations 
with the parallel Vlasov code LOK l 26 i 32 of a LW driven 
with an external potential with a smooth transverse pro- 
file. Our results are in qualitative agreement with Sec. 
IV of Ref. [2^. This work considered a free LW excited 
by a driver of finite duration, while we consider a driver 
that remains on. 

We present a unified calculation of collisional detrap- 
ping due to electron-ion and electron-electron collisions, 
including both pitch-angle scattering and parallel slowing 
down and diffusion. This relies on the observation that 
the diffusion rate increases rapidly with Fourier mode 
number in velocity, so that only the electrons in the fun- 
damental mode remain trapped after a short time. The 
collisional detrapping rate scales as 1/SN, since the trap- 
ping width in velocity increases with wave amplitude. We 
discuss two ways to compare the relative importance of 
detrapping by side loss and collisions, which is compli- 
cated by their different scaling with 5N. 

Our calculations are applied to ICF plasma condi- 
tions, particularly LWs driven by stimulated Raman 
backscattcr (SRBS) on the NIF. Transverse side loss out 
of laser speckles in a phase-plate-smoothed beam is gen- 
erally a more effective detrapping process than collisions. 
The threshold 5N for trapping to overcome side loss 
decreases with density and increases with temperature, 
while the collisional threshold decreases with density and 
slightly increases with temperature. For conditions typ- 
ical of backscatter on NIF ignition experiments, namely 
n e = 0.1n cr and T e =2 keV, a reflectivity of (5.28 x 10 13 
W cm~ 2 //o) 2 produces linear Langmuir waves above the 
side loss threshold. Such values are likely to occur in in- 
tense speckles of phase-plate-smoothed NIF beams. We 
also show that smoothing by spectral dispersion (SSD) 33 
is ineffective at detrapping. 

The paper is organized as follows. Section HTl provides 
some general considerations on our detrapping analysis. 
We present in Sec. IIIII convective loss calculations for 
both longitudinal (end) and transverse (side) loss. Sec- 
tion |IV] contains Vlasov simulations with the LOKI code 
which study the competition of trapping and side loss. 
Detrapping by Coulomb collisions is treated in Sec. [V] 
Our results are applied to SRBS in underdense ICF con- 
ditions in Sec. EU We conclude in Sec. IVII1 and present 
some mathematical details of our collisional derivation in 
the Appendix. 



II. GENERAL CONSIDERATIONS 

Consider the group of resonant electrons in a LW field, 
attempting to undergo bounce orbits. There is first a 
time-dependent condition for trapping to be significant, 
even in the absence of any detrapping process. For in- 
stance, if a LW is suddenly excited in a Maxwellian 
plasma, electrons execute bounce orbits according to 
what we call the dynamic bounce number 

The time dependence of tb allows for a slowly- varying 
wave amplitude Sn(t). We use the subscript x since the 
longitudinal coordinate x is the direction of wave prop- 
agation. Vlasov simulations presented in Sec. |IV] show 
that trapping significantly affects the dielectric response 
when Nb,x ~ 0.5. That is, it takes a finite time for the 
dynamics of resonant electrons to reflect trapping. The 
early works of Morales and O'Xci! 2, ' indicate such behav- 
ior, where the damping rate and frequency shift evolve 
over several bounce periods until approaching steady val- 
ues as the system reaches a Bernstein-Greene-Kruskal 
(BGK) state!. 

To estimate the threshold for trapping to overcome a 
detrapping process, we assume the wave has been present 
long enough that Nb,x is large, and we can treat all res- 
onant electrons as trapped. For coherent trapping effects 
to occur, an appreciable fraction of electrons in the trap- 
ping region at some time must remain there for about 
a bounce period before being detrapped. We are inter- 
ested in the number of electrons in the trapping region, 
and how long they stay there. 

We define the "trapping region" to extend from u = 
u p ±u tr /2 where u± r — A{k\Be)~ 1 SN 1 / 2 is the full width 
of the phase-space trapping island and Xdc = VTe/uJ pe 
with VTe = (Te/me) 1 / 2 . Throughout this paper, we use 

ux = vx/vTe (2) 

to denote the scaled velocity vx for various subscripts X. 
Let Nt T (t) denote the fraction of electrons in the trapping 
region at the initial time t = 0, that continuously remain 
so to some later time t (note N tT (t = 0) = 1). At t = 
we take the electron distribution to be Maxwellian. The 
fact that only some electrons in the trapping region lie 
within the separatrix (depending on their initial phase 
kx) is not relevant here, since all the detrapping processes 
considered here are insensitive to the electron's phase in 
the wave. That is, the rate at which electrons leave the 
trapping region is independent of kx. 

The detrapping rate is defined by assuming exponen- 
tial decay at the rate Vd = l/td for the trapped fraction: 
iVtr = e~ Udt . We allow for several independent detrap- 
ping processes to occur simultaneously, in that the overall 
detrapping rate va.o is the sum of the rates — 
for each ith process considered separately. Since a de- 
trapping process generally does not strictly follow expo- 
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nential decay, we choose a critical fraction 7V t * r , which ob- 
tains for a critical time t — t* , and let v d — m(l/iV* r )/i*. 
This Vd is independent of iV t * for exponential decay. We 
shall choose N* T = 1/2 in what follows. Given the ap- 
proximate nature of our calculation, further refinement 
of v c i has little value. 

We sometimes make use of a ID kinetic equation with 
a Bhatnagar-Gross-Krook relaxation (or simply a Krook) 
operator- 
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[d t + vd x - (e/m e )Ed v ] / = v K ■ (n/ /n - /). (3) 

The Krook operator relaxes the electron distribution 
function / to an equilibrium / , and locally conserves 
number density n = J dv f. Such a model is sometimes 
used to mimic detrapping by transverse convective loss 
(a higher space-dimension effect) or Coulomb collisions 
(a higher velocity-dimension effect), such as in Ref. l35l . 
Any perturbation from fo decays exponentially at the 
rate vk, so Vd = Vr for such an operator. This is espe- 
cially useful for a detrapping process which has Vd inde- 
pendent of wave amplitude; this is the case for convective 
loss but not for collisions (as shown below). SRS simu- 
lations with a ID Vlasov code and Krook operator, and 
its suppression of kinetic inflation, are presented in Ref. 

We take the bounce period of all trapped electrons to 
be tb, the result for deeply-trapped electrons. The ac- 
tual period slowly increases to infinity for electrons near 
the separatrix. We then define the bounce number for 
process i as 



N BA = — 



TB 



5N 



SN 



(4) 



We have expressed Nb,i as a ratio of the LW amplitude 
to a "threshold" amplitude 5Ni to some power pi. Re- 
call that trapping effects like the Landau damping re- 
duction develop continuously with SN, so the thresh- 
old for trapping nonlinearity is not a hard one. Be- 
sides the 8N~ X / 2 dependence of tb, td,i also depends 
on SN in a process-dependent way. The overall detrap- 
ping rate v^o = J2i v d,i, gives an overall bounce number 
N B: o = t d ,o/TB which can be written N~^ Q = £\ N B \. 
We also define an overall threshold amplitude S No such 
that Nb o[$N = SNo) = 1; it is not generally true that 
SN Q V-VY,. 



III. CONVECTIVE LOSS: THEORY 
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FIG. 1. (Color online.) Trapped electron fraction -/V tr , a i due 
to transverse side loss for a 2D plane (blue) and 3D cylinder 
(red) region. The dashed curves are the appropriate early and 
late time limits. See Eqs. ([5]) through 



For v tr v p we can treat all the electrons in the trap- 
ping region as moving forward at v p . Thus for end loss 
N. 



tr,e] 



1 - v p t/L\\. To find i/d,el) we take N* x cl = 1/2, 
L\\/v p and = K c \v p /L\\ with 



which gives t* x 
A' i = ln2. 

For transverse side loss, consider a cylindrical 
wavepacket of transverse diameter L± and infinite lon- 
gitudinal length. In N total spatial dimensions, the 
cylinder has an N — 1 dimensional cross-section. Elec- 
trons with a Maxwellian distribution are transiting 
the cylinder, with unnormalized distributions f± = 
u^ 2 cxp[— tij_/2] where u±_ is the transverse speed and 
f±du± is the number of electrons per du±. The average 
u± = ([tt/2] 1 /2 ; [8/tt] 1 / 2 ) = (1.25, 1.60) for N = (2, 3), in- 
dicating the detrapping rate should be faster in 3D than 
2D. 

In any dimensionality N , is a universal function 

of t = tvxe/L±. In 2D we find 



l/i 



A t 2 r D sl = (2n)-^ / du e"" '\l - \u\t) (5) 

J-i/i 

= erf [l/IV2l + (2/^) 1/2 t(e- 1/2p - 1). (6) 



In a LW of finite spatial extent, electrons remain in the 
trapping region only until they transit the wave. This 
detrapping manifests itself by longitudinal loss out of 
the ends of the wavepacket (the x direction for our field 
representation cos(fcir — tot)), as well as transverse loss 
out the sides. End loss is simply found by considering a 
wavepacket of length Lu and infinite transverse extent. 



In 3D we obtain 



A 3D , 

Jv tr,sl 



i/i 



du u e~" 2/2 



1 (arcsin[ut| + ut[l - (ui) 2 } 1/2 ) 

7T 



(7) 
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The limiting forms are 



N, 



2D 



(i»l) 



/v 2D , 



^t 3 r D sl(^«l) 
^t 3 r D sl(^»l) 



1 - [2/^/% 

l-[Z/-K] i m, 

1/8?. 



(8) 
(9) 

(10) 
(11) 



In both limits the decrease is more rapid in 3D than in 
2D. Figure [1] displays the various formulas for N tIlS \(t). 

The resulting detrapping rates, based on iVtr,si = 1/2, 
are z/ d>s i = K sl v T e/L± with K sl = (1.02, 2.08) in '(2D,3D). 
As expected, the 3D detrapping rate is faster. The 3D de- 
trapping rate exceeds the 2D one by a larger factor than 
the average transverse speed because the faster electrons 
leave first, and the relative surplus of electrons in 3D 
over 2D (proportional to u±) increases with transverse 
speed. A wavepacket with asymmetric (e.g. elliptical) 
cross-section should have a rate between the 2D and 3D 
result with L± taken as the shortest transverse length. In 
a laser beam smoothed with phase plates, elliptical speck- 
les can be produced by certain polarization-smoothing 
schemes or a non-spherical lens; Langmuir waves driven 
by SRS in such speckles would also acquire an elliptical 
cross-section. 

Comparing the end loss and side loss rates gives 



(12) 



For the LW to not experience strong Landau damp- 
ing, we have v p > vtc- Of course, L±/L\\ depends on 
the physical situation (laser speckles are discussed in 
Sec. IVI|) . The bounce number for side loss is Nb, s \ = 
[SN/SNsi] 1 / 2 , with exponent p s i — 1/2 and threshold 
amplitude SN S \ = [2irK s \A[) e / L±] 2 . In practical units 
5N S ] = (9.44 x 10 18 /n e , cc )T^ kv / L± jlMm where n e . cc is in 
cm~ 3 , Te.fcv is m keV, and L±^ m is in /im. The end loss 
result is completely analogous: Nb, c \ = [SN/SNgi] 1 ^ 2 , 
8N cl =[27rK cl u p X De /L^ 2 . 



IV. VLASOV SIMULATIONS OF CONVECTIVE SIDE 
LOSS 

In this section, we quantify the competition between 
convective side loss and electron trapping in a driven 
Langmuir wave. We use the parallel, 2D-2V Eulerian 
Vlasov code loki 32 . This code employs a finite- volume 
method which discretely conserves particle number. The 
discretization uses a fourth-order accurate approximation 
for well-resolved features, and smoothly transitions to a 
third-order upwind method as the size of solution features 
approach the grid scale. This construction enables accu- 
rate long-time integration by minimizing numerical dis- 
sipation, while retaining robustness for nonlinearly gen- 
erated high frequencies. As a result, the method is not 
strictly monotone- or positivity-preserving, nor does it 
eliminate the so-called recurrence problem. Recurrence 



occurs at a recurrence time of t rec = A/Au when further 
linear evolution of a sinusoidal perturbation cannot be 
represented on a given grid. 

Our simulations are electrostatic with fixed ions and 
an external driver electric field E^x: 



E d = E A(t)h(y) cos(k x - uj t). 



(13) 



There is no y component to the driver field, which would 
be needed if the driver were derived from a scalar po- 
tential. We shall only discuss the x component of the 
field. The temporal envelope A(t) ramps up from zero to 
unity over a time 50/cj pe and then stays constant. The 
transverse profile h(y) is 



2ny 



|(1 + cosk iy ), \y\ < (14) 



otherwise. 



(15) 



fci = Ait/Ly. All runs were performed over one driver 
wavelength in x, with periodic boundaries for fields and 
particles, and N x = 32 zones in x (note that runs 
performed with N x = 64 gave qualitatively similar re- 
sults). Runs with a transverse envelope h(y) had peri- 
odic boundaries for fields and particles at |y| = L y /2, 
and a Krook operator with a y-dependent damping rate 
to rapidly relax the distribution to the initial Maxwellian 
near the boundaries. The runs were thus effectively fi- 
nite in y. We used N y = 11 to 45 zones in y, with more 
used for larger L y and to check convergence. The v x and 
v y grids both extended to ±7«Te- N vy — 32 zones in v y 
were used throughout. N vx is set by two requirements: 
the trapping region must be adequately resolved, and re- 
currence phenomena must not be significant. We found 
Av x ~ 0.1t>i r was sufficient to give converged results. 
loki's advection scheme is designed to mitigate aliasing 
problems, and we only saw modest effects related to it 
when comparing runs with different N vx . 

We first present ID runs with h(y) — > 1, which are 
detailed in Table HI The total electric field magnitude Et 
(drive plus internal) from linear theory is 



Et 



1 



1 + X 



[(1 + Rc X ) 2 + Im X 2 ] 



-1/2 



(16) 



The electron susceptibility x appropriate for the Krook 
operator of Eq. 



x(w,fc) 



z'(C) 

2{k\ De ) 



vk 



kVTeV2 



z(0 



(17) 



where £ = ui/kvTeV2 and Z is the plasma dispersion 
function^. For analyzing our LOKI runs we set vk = 0, 
but use vk below for collisional damping of LWs. \ 
is evaluated at the driver ko and ujq. We chose ujq 
to give nearly the maximum Et for a given fco. For 
ko^De < 0.53, a resonant ujq exists where 1 + Rex = 0; 
the maximum Et then occurs close to this point. No res- 
onance exists for koXne > 0.53, which is called the loss 
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of resonance 3 -^. Of course some u>o still maximizes Et 
in this regime. The non-resonant case differs from the 
resonant one, in that reducing Imx and Landau damp- 
ing, e.g. by flattening the distribution at the phase ve- 
locity by electron trapping or some other means, does 
not lead to a large enhancement in the Langmuir wave 
response to an external drive. The term 1 + Rex in Eq. 
(1161) keeps Et hnite even if Im^ = 0. For the parame- 
ters of the run ID. 7a, we find Et/E = 1.80 for the full, 
complex x-i while setting Im% = slightly increases it to 
E T /E = 2.01. 

Similar logic applies to kinetic inflation of stimulated 
Raman scattering. Electron trapping and the resultant 
Landau damping reduction can greatly increase the scat- 
tering at a resonant wavelength. However, scattering 
at a non-resonant wavelength is not subject to infla- 
tion, and can even decrease, due to reducing Im%. Non- 
resonant SRS can occur in a situation seeded away from 
resonance^, or if the plasma conditions are such that no 
resonance exists for any scattered wavelength, namely 
high T e and low n e . 

Figure [2] presents the results of our ID runs. The left 
panel shows the time evolution of the total electric field 
amplitude, normalized to the linear value from Eq. (1161) . 



Early in time (ui pe t = 100 — 200) the linear response is 
achieved, which validates the linear dispersion and prop- 
erties of loki when using the chosen grid resolution. As 
time progresses the response increases due to the damp- 
ing reduction, and then oscillates due to the interplay 
of the frequency shift and the fixed driver. Similar be- 
havior was seen in Ref. [28|. We plot the results vs. the 
dynamic bounce number Nb, x from Eq. (fTJ), using the 
time-dependent total field, in the center and right pan- 
els. Nb, x is thus a trapping-based re-scaling of time. 
The other runs from Table [T] are included as well. The 
driver strength Eq was chosen in runs ID . 35b, ID . 5a, 
and ID. 7a to give similar bounce periods. In all cases, 
the linear response is achieved after a transient period 
related to driver turn-on, until Nb, x ~ 0.5. After this 
point the response increases, until the frequency shift 
develops at Nb,x ~ 1- As koXoe increases, the enhance- 
ment above linear response decreases. This is likely due 
to the rapid increase of the frequency shift with kXr> e , as 
shown by most theoretical calculations, e.g. Ref. [j. For 
ko^-De = 0.7, there is a slight enhancement to 1.3x the 
linear response, followed by a dip to about 0.7x and sub- 
sequent oscillation about unity. This lack of significant 
trapping nonlinearity agrees with the above discussion of 
the non-resonant regime. 

Following the side loss analysis of Sec. IIII1 the 2D side 
loss rate is z^si = 1.02vTe/ L±. We take L± — L y /4, 
the full- width at half- max of h(y). The side loss bounce 
number is then 



25.6\ De 



-8N 1 ' 2 . 



(18) 



Gauss's law relates 6N to the total electric field: SN — 
fcoAoe ■ Et where E = Ee/m e VTeU P e is the normalized 



TABLE I. ID loki runs with no transverse driver profile. 

Eq — Eoe/m e VT e U3pe- 



Run fcdA.De Ulo/lOpe Et/Eq 



E 



TBUpe plot curve 



ID. 35a 0.35 
ID. 35b 
ID. 35c 

ID. 5a 0.5 

ID . 5b 

ID. 7a 0.7 



1.22 



11.9 



1.44 3.22 



1.79 



1.25X10 -3 
5xl0" 5 
2xl0~ 4 
1.3xl0~ 4 
5.2xl0" 4 
1.7xl0 -4 



871 red dash 

436 solid black 

218 blue dash 

434 black dash 

217 green dash 

430 black dot 



TABLE II. 2D loki runs with transverse driver pro- 



file. All 


runs have 


ko^De — 0.35, loo/uj 


pe = 1.22, and 


E e/m e v T 


eUipe = 2 x 10 , the same as run 


ID. 35c. 


Run 


Ly/\De 


Nb,a = W884A De 


plot curve 


2D100 


100 


0.113 


solid red 


2D200 


200 


0.226 


solid dark blue 


2D400 


400 


0.452 


solid green 


2D800 


800 


0.905 


solid magenta 


2D1200 


1200 


1.357 


solid blue 



electric field. Using the linear response from Eq. (|16[) . we 
obtain the linear estimate 



N, 



B,sl 



25.6\ De 



ko^De 



1 + X 



1/2 



E, 



1/2 



(19) 



The 2D loki runs are listed in Table HU All runs used 
k \ De = 0.35, Lo/uj p e = 1.22, and E = 2 x 10~ 4 , the 
same as run ID . 35c. For these values, our linear estimate 
becomes Nb, s \ — L y /884A,De- 

The total electric field magnitude, E T {y = 0), is plot- 
ted vs. the dynamics bounce number Nb x (also found 
using the y = field magnitude) for our 2D runs in 
Fig. [3] The black curve is the analogous ID run. For 
Nb,x ^$ 4 there is a continuous increase in the response 
with profile width L y . This allows us to quantify trap- 
ping nonlinearity vs. L y , which we do in Fig. 3J The 
abscissa in that figure is the side loss bounce number, 
Nb, si, computed with linear response as in Eq. (|19p . The 
ordinate is the field enhancement due to trapping, scaled 
to the same quantity for the ID run. This is shown at 
times corresponding to Nb, x — 1 and 2, which are both 
early enough that the amplitudes have been mostly in- 
creasing, with little oscillation due to the frequency shift. 
The two curves agree well, and demonstrate the contin- 
uous development of trapping effects with wide profiles. 
Slightly more than half the ID trapping effect obtains 
for Nb. s \ — 1, which vindicates our Nb ~ 1 approximate 
threshold for trapping. 

The plasma response to a driver with transverse pro- 
file h(y) differs from the ID case. This can be seen in 
the ordinate of Fig. 2] falling below zero for the small- 
est L y — 100. There have been several linear calcula- 
tions of transit-time damping in LWs of finite extent, 
mostly by integration along particle orbit o 38 ' 39 . Ref. l38l 
showed that, for a potential with a step- function profile 
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FIG. 2. (Color online.) Left: Et = magnitude of total electric field (x component), scaled to linear response, vs. time for ID 
LOKI run ID. 35b. Center: Et vs. dynamic bounce number Nb, x from Eq. (JXJ) using Exit) from ID loki runs in Table [Q (see 
table for curve meanings). Right: center panel for expanded domain. 
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N B ,x 

FIG. 3. (Color online.) Magnitude of E T {y = 0) for 2D loki 
runs with transverse driver profile h(y) of various L y . Run 
parameters and curve meanings are given in Table HU Black 
solid curve is ID run ID. 35c. 



in space, the transit-time damping exceeds that for an in- 
finite plane- wave for uj/kvTe 3> 1, while for Lo/kvxe ~ 1 
it can be less. We adopt the alternative approach of 
writing the response as a superposition of responses to 
the Fourier modes comprising the drive. This is particu- 
larly convenient for our h(y), which is composed of only 
two Fourier modes. For simplicity we present the result 
for h(y) periodically repeated, instead of the actual LOKI 
profile with compact support over \y\ < L y /A. The com- 
pact case would lead to a continuous Fourier transform 
rather than discrete series, and introduce a line width 



[Ex-Exjin] / [Ex 1D -E x ,lin] 
-| o J i 1 1 1 1 1 1 1 1 J 1 1 1 J 1 1 1 J 1 1 1 J 1 1 1 J 1 1 1 J_ 




-| i i i i i i i i i I i i i i i i i i i I i i i i i i i i i p 

0.0 0.5 1.0 1.5 

Nb,sI 



FIG. 4. (Color online.) Departure from linear response for 
2D LOKI runs, scaled to the same quantity from the ID run 
ID. 35c. The (black, red) curves are taken at a time when the 
dynamic bounce number Nb,x is (1,2) respectively. Nb, s i is 
found from linear response, using Eq. Q19fl. 



around the dominant modes. This does not change the 
qualitative result. Unlike Ref. l38l . our compact profile 
h(y) is not a step function but smooth, with h and h' 
continuous at all points (although h" is not). 
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The drive Ed, made periodic in y, is 



E d 



1 4. _Jkiy , _ f ,—ikiy 

+ 2 + 2 



+ c.c. (20) 



A standard kinetic calculation, accounting for the fact 
that the LOKI driver field has no y component and thus 
does not come from a potential, gives the total field at 
y = 0: 

E T (x,t,y = 0) = E \R\ cos(fc x - uj a t + a), (21) 
1 , l + (l + (fco/fc 1 ) 2 )- 1 x+ 



2R = 



1 



Xo 



1 + X+ 



(22) 



a is a real phase. \ is the usual collisionless susceptibility, 
which depends only on oj and k = \k\. xo — x(ko,u>o) and 
X+ = x{k+,uo) with k + = (fcg + k 2 ) 1 / 2 . For k\ — 0, we 
recover the ID result Eq. (fT6|) . Physically, the higher- fc 
modes induced by the transverse profile are more Lan- 
dau damped (as well as being slightly off resonance for 
the fixed ujq), which reduces the response. For the pa- 
rameters of Table HI we find = (0.801,0.948) 
for L y = (100,200) where \R\ 1D = 11.9 is the value for 
L y —> co. We thus find a slight decrease in the lin- 
ear response for our sharpest profile (L y = 100), and 
an insignificant change for wider ones. This is borne 
out by Fig. [3] The red curve for L y = 100 shows no 
signs of trapping, and reaches a steady level slightly 
more than 0.8 times the ID linear value. The blue curve 
{L y = 200) shows a slight trapping enhancement, and 
reaches a steady level slight above 1.2x linear after about 
2 bounce periods. 



V. COULOMB COLLISIONS 

Collisions remove electrons from the trapping re- 
gion via pitch-angle scattering (from electron-ion and 
electron-electron collisions) as well as parallel drag and 
diffusion (from only electron-electron collisions since 
mi/m e 3> 1). We adopt a Fokker-Planck collision op- 
erator: 

d t f = + Z cS )u~ 3 dp [(1 - M 2 )9 M /] 

+2is u- 2 d u {f + u- 1 d u f). (23) 

u = cos 9 where 9 is the pitch angle between u and the u x 
direction, and u = \u\. Vq vs, a. thermal electron-electron 
collision rate: 



_ W pe lnA e 



87r7Vr 



(24) 



Nn e = n eX 3 De and In A ee = 24 — \n(nl^ 2 /T e ) (n e in cm 3 , 
T e in eV) is the electron-electron Coulomb logarithm ap- 
propriate for T e > 10 eV (Ref.ffl p. 34). The effective 
charge state is 



^off 



\- fjZj lnA e 
^ Z lnA e 



(25) 



where nj = J^i n i i s the total ion density, Z — J^i Zifi 
with fi = rii/nj; J^. fi — 1, and In A ei is the electron- ion 
Coulomb logarithm^. 

In section IVII we apply our results to Langmuir waves 
generated by Raman scattering in underdense ICF plas- 
mas, which are typically low-Z. For instance, NIF igni- 
tion hohlraum designs currently use an He gas fill (with 
H/He mixtures contemplated), and plastic ablators (57% 
H, 42% C atomic fractions). This gives Z c g = 5.08 when 
fully- ionized and lnA ei = lnA ee . Be and diamond abla- 
tors are also being considered. For illustration, we take 
Z c ff = 1 as the lowest reasonable value (fully- ionized H), 
and use Z g = 4 (fully-ionized Be) to represent an abla- 
tor plasma. 

It is useful to define a unitless time t (different from 
the side loss t used above), which demonstrates some of 
the basic collisional scaling: 



In 1 

tt 2 is 



(kXi 



tt (k\ De ) 5 lnA e 
128 (w/w pe ) 3 N De 



(26) 
(27) 



Our calculation of the trapped fraction is detailed in the 
Appendix. The key observation is that electrons not in 
the fundamental Fourier mode in u x diffuse very quickly 
out of the trapping region, so it suffices to consider just 
the number in the fundamental mode. At t = 0, this is 
81%. The upshot is that A^ r ,cj the fraction of initially 
trapped particles remaining in the fundamental mode af- 
ter time t, is 

/•OO 

iV t r,c(*, Z cS , u p ) = 0.81 / du ± u ± exp [-u 2 x /2 - Dt\ . 
Jo 

(28) 

D(u±, Up, Z e g) is given in Eq. (|A12[) . 

Eq. |28|) is an implicit, integral equation for t as a func- 
tion of Z e g, Up, and iVt r . c . We find the "exact" solution 
by performing the integral numerically, and interpolating 
t for a desired N tT , c - We derive an approximate solution, 
valid for u p 3> 1, for t in the Appendix. The result is 



t ~ tn + i\u r . 



(29) 



to and t\ are both positive, and the decay time decreases 
with increasing u p . Figure [5] plots A triC (£) for several 
Up and Z c ff, using the exact results (solid curves) and 
the approximate form for u p — > oo of Eq. (|A18j) (dashed 
curves). Few electrons remain trapped at t = 1. The 
approximate forms are quite good, even though u p is not 
that large. 

Figure [S] displays the relative error e = 1 — t appi: /t ex 
between t for iVtr.c = 1/2 computed two ways. The exact 
t ex is found numerically, and i a ppr is from Eq. (|29[) . with 
Eq. (|A19|1 for t and Eq. (|A24p for ii. The agreement is 
excellent, within 1% for most of parameter space. 

The collisional detrapping rate v^.c is 



v c ln(l/jV tr ) 
6N t 



(30) 
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J= 0.1 - 




100.0 



t = v c t / 5N 



Note that 2/<j,c ~ SN^ 1 since u tr ~ SN^ 1 / 2 . The reason is 
that the larger the wave amplitude the wider the trapper 
region extends in velocity, and collisions take longer to 
remove the electron velocity from this region. Recall that 
u^c depends slightly on the choice of N tr due to the non- 
exponential decay of N tT with t; as with convective loss 
we choose N tI — 1/2. 

The collisional bounce number is 



N B ,c = 



SN 



-i 3/2 



5N C 



6N r = 



2nv. 



ln(l/JV tr ) 
i 



2/3 



(31) 

The amplitude exponent for collisions is p c = 3/2, unlike 
the convective loss value of 1/2. This stems from the 
fact that Vd for collisions is amplitude-dependent while 
for convective loss it is not. We now construct the overall 
bounce number Nb,o for convective side loss and colli- 
sions, as outlined above. Assuming that separate de- 
trapping processes are independent, and their detrapping 
rates add, yields 



N 



B,0 



= N 



B.sl 



\SN sl - 


1/2 


~6N C ~ 


3/2 


SN 


+ 


SN 





(32) 



We define an overall threshold amplitude SNq such that 
Nb,o[5N — 5 Nq\ — 1. Eq. (|3"2")l gives a cubic equation 





Up = V P / VTe 

FIG. 6. (Color online.) Relative error e = 1 — t apP r/tex be- 
tween the exact and approximate results from Eqs. (|28I29[) . 
respectively, for collisional t; see text for details. The green 
curves are e = [-0.01, -0.01, 0.05]. 



FIG. 5. (Color online.) Trapped fraction due to collisions 
iVtr.c vs. unitless time t defined in Eq. ((26}. Solid curves are 
exact results from Eq. pg]l. (u v , Z cB ) = (2,1), (4,1), and (4,4) for a : 
for black, red, and blue, respectively. Dashed red and blue 
curves are approximate results for u p —> co from Eq. (TM8|) . 
which depend only on Z B g and not u v . Green dashed curve 
is 0.81 exp[— 2t\, the approximate form neglecting the term 
proportional to t in the denominator. N t r, c {t = 0) = 0.81 and 
not unity due to electrons not initially in the fundamental u x 
mode. 



5N /2 : 



5Nll\. 



5N? /2 = 0. 



(33) 



There are two ways to compare the relative impor- 
tance of side loss and collisions. One is: for which pro- 
cess must the wave amplitude SN be larger for trapping 
to be significant (Nb = 1)? The other is: for a given 
SN, which process will detrap more effectively? The two 
views are not equivalent, due to the different dependence 
of the side loss and collisional detrapping rate on 6N. 
The first amounts to comparing the thresholds 5N S \ and 
5N C , which can be computed just from plasma and wave 
properties without knowing SN. The ratio of detrapping 
rates can be written in terms of a critical amplitude SN cr : 



SNcr 

SN 



SN rr = 



In 2 



tK s \ W pe Xoe 



(34) 



VI. PARAMETER STUDY FOR ICF UNDERDENSE 
PLASMAS 

We now apply our analysis to ICF conditions where 
stimulated Raman scattering (SRS) can occur, namely 
the underdense coronal plasma. SRS is a parametric 
three-wave process where a pump light wave such as a 
laser (we which label mode 0) decays to a scattered light 
wave (mode 1) and a Langmuir wave (mode 2). We 
restrict ourselves to exact backscatter (SRBS; k\ anti- 
parallel to ko), as this generates the largest (small- 
est v P 2/vtc) and thus makes trapping effects more im- 
portant (small transverse components to k 2 have little 
effect on the phase velocity). Both measurements and 
simulations with the paraxial-envelope propagation code 
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pF3D^i have shown backscatter to be the dominant di- 
rection for SRS. With fcj = kiZ, the phase-matching 
conditions are loq = uj\ + ^2 and fco = k\ + &2 with 
fci < 0. We employ the (cold) light-wave dispersion re- 
lation w? = (cfci) 2 + Wp e for modes i =0 and 1, and use 
the vacuum wavelength Ai = 2nc/uji. Frequency match- 
ing thus requires n e < n cr /A, with n cr ^ = (eom e /e 2 )u; 2 
the critical density for mode i, and n cr — n cr _Q. For spe- 
cific examples we choose Ao= 351 nm, appropriate for 
frequency-tripled UV light currently in use on NIF. Spe- 
cific plasma conditions thought to be typical for SRBS 
on NIF ignition targets, during early to mid peak laser 
power, are n e /n cr = 0.1 and T e = 2 keV (Ai ~ 550 
nmp. The scattered wavelength continuously increases 
during a NIF experiment, consistent with the hohlraum 
filling to higher density. 

An important case for this paper is LW's driven by 
SRBS in the speckles of a phase-plate-smoothed laser 
beami^. For a laser wavelength Ao and square RPP with 
optics F-number F, the intense speckles have L± ps FAo 
and L|| ps 5F 2 Ao (see Ref. A speckled beam is 

not the only situation where SRS can occur; for in- 
stance, there has been recent interest in re-amplification 
of backscatter by crossing laser beams 4 ^ and backward 
Raman amplifiers 4 ^. However, for a single laser beam, ex- 
periments at Omega and pF3D simulations show speckle 
physics, and its modification by beam smoothing, must 
be accounted for to accurately model SR S 47 ' 48 . Experi- 
ments have also verified the increase in backscatter with 
increased gain per speckle length, by changing the laser 
aperture and thus the effective F— . We therefore focus 
on speckles. On NIF, four laser beams, each smoothed 
by a phase plate and with an overall F = 22 square aper- 
ture, are grouped into a "quad" which yields an effective 
square aperture of F pa 8. We thus use F — 8 for illus- 
tration. As the beams of a quad propagate through a 
target, they can separate from one another, refract, and 
undergo other effects that change the shape of their ef- 
fective aperture and speckle pattern. We do not pursue 
this further here, but it should be born in mind when 
applying our analysis. Also the ratio L±/Ln — 1/5F is 
so small that ^d,el/^d,sl ~ {K e \/ K s \)u p /bF = 0.0083u p 
(3D) is small for essentially all speckles of interest. Thus 
side loss is a more potent detrapping mechanism than 
end loss, in speckles. 

To quantify detrapping rates, we consider the thresh- 
old amplitudes 8N S \ and 8N C . Unlike SN B i, SN C depends 
on L02 and &2 of the Langmuir wave. For a given set of 
plasma conditions, the choice of (0^2,^2) is not unique 
but depends on the application. For SRS developing lo- 
cally, one can choose the LW corresponding to the largest 
growth rate for those conditions. Another approach is to 
consider a single scattered-light frequency as it propa- 
gates through a target. We consider only k variations 
induced by spatial profiles and not oj variations due to 
temporal plasma evolution^ (which is mostly relevant to 
stimulated Brillouin scattering) . In this case, the match- 
ing conditions given the local plasma properties dictate 
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FIG. 7. (Color online.) Wavelength Ai in nm for SRBS light 
(increments of 50 nm), for a pump wavelength Ao= 351 nm. 
Black solid: Ai phase-matched with a natural Langmuir wave, 
satisfying the dispersion relation Eq. (|35p with vk = 0. Red 
dash: Ai for the maximum local SRBS spatial gain rate. 
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FIG. 8. (Color online.) Langmuir wave feA^e phase-matched 
for SRBS with different Ai choices, for a pump wavelength 
Ao= 351 nm. Black solid: Ai phase-matched with a natural 
Langmuir wave, as in black solid curve of Fig. [7] Red dash: 
Ai=550 nm (cj's fixed, k's vary). Blue dot: Ai=550 nm black- 
solid contour from Fig. [7] 
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FIG. 9. Threshold Langmuir wave amplitude for side loss, 
5N„i, for L ± =8x351 ran. 



how k 2 varies. 

Figure [7] presents the local Ai for SRBS computed in 
two ways. The black curves are found by phase-matching 
with a "natural" LW, by which we mean uj 2 = Rc[oj2 C ] 
where complex uj 2c satisfies 



1 + X[ fc 2r,W 2c ] = 



(35) 



with real k 2r = k§ — k\. To find lu 2c , we set vk = and 
recover the usual collisionless X- We use vk ^ below 
as a simple way to include collisional LW damping when 
Landau damping is negligible. The red curves in Fig. [7] 
are the Ai which maximizes the local spatial SRBS gain 
rate in the strong damping limit^i: 



5 2 lnii(Ai,z) 



2irr e I 



m e c 2 uoko 



\ki\ 



Tm- 



1 



(36) 



We use the collisionless x with vk = 0. The first bracket 
is independent of Ai, while the second bracket does. The 
two results for Ai in Fig. [Jj are very close except for high- 
fc2A.De LW's (low n e , high T e ), where Landau damping 
and its variation with Ai is significant. We choose for con- 
venience to use Ai matched to a natural LW below. We 
display in Fig. [5] the k 2 X De corresponding to two choices 
of Ai. The black curves use the Ai phase-matched to a 
natural LW (the black curves in Fig. [7J, while the red 
curves are for a constant Ai= 550 ran. 

The side loss threshold <5iV s i is shown in Fig. |H1 for 
L± = FAq and F = 8. It simply represents the varia- 
tion in Adc, and is independent of (uj 2l k 2 ). Figure [TOl 
depicts the collisional threshold SN C for Z c g= 1 and 4. 
The decrease of SN C with electron density is mainly due 
to the decrease of the (kXne) 5 factor in v c (see Eq. 
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FIG. 10. (Color online.) Threshold Langmuir wave ampli- 
tude for collisions, SN C . LW 0J2 and fe are for phase-matched 
SRBS with a natural LW (black solid curve of Fig. [7J, and 
pump wavelength Aq= 351 nm. Black solid: Z e g=l. Red 
dash: Z e ij=4. From right to left, curves are for SN C = 
[1,2,5,10,30] x 10~ 4 . 

which in turn is due to the 1 /Vp fall in the Coulomb cross- 
section (see Eq. (|24|) L The ratio 5N c /SN s i is displayed 
in Fig. Ill( which indicates collisions have a minor effect 
except for low n e and low T e ; this relative importance 
depends strongly on the transverse length L± chosen for 
side loss. Figure [12] plots the critical amplitude 8N cr 
from Eq. (|34|1 . For SN > SN cr , the side loss detrapping 
rate exceeds the collisional rate. SN cr is larger at smaller 
n e , indicating collisional detrapping is more relevant. 2D 
particle- in-cell simulations with the VPIC code of Raman 
amplifier experiment o 45 ' 52 found that collisions mattered 
for low-intensity seed light waves in a low-density plasma 
(n e /n cr — 0.01). 

The reflectivities which correspond to trapping nonlin- 
earity can also be estimated. We assume the LW's are in 
the strong damping limit^, and write 



SN = 



1 (k 2 \ D e) 2 V V 1 



2 1 



(37) 



Vi = eEi/m e LUi is the oscillation velocity for light wave i\ 
in practical quantities we have (Vi/c) 2 = IiXf /ijiPo, with 
P = 2n 2 (e /e 2 )m 2 e c 5 = 1.37 x 10 18 W cm~ 2 • ^m 2 . 77; = 
[1 — rie/ncr^] 1 / 2 reflects the decrease in group velocity. 
With reflectivity R — Ii/Iq, we find 



Ir.r 



(k 2 Xi 



,) 2 m e c 2 A Aj 



-1/2 



(38) 
(39) 
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FIG. 11. Ratio 5N c /5N Bl from Figs. [5] and [TO] for Z cS =4. 
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including the Krook operator vk to damp the LW 
when k2Xoe is small and Landau damping is ineffective. 
For this purpose we choose vk to be the collisional, un- 
magnetized frictional drag rate in the electron momen- 
tum equatio n 53 ' 54 . This rate is appropriate for the drag 
on the bulk sloshing motion of the electrons in the LW 
electric field, and not collisions of resonant electrons with 
u « Dp. I cr minimizes near the lower-right corner near 
the 2.5 x 10 14 curve, where Landau and collisional damp- 
ing are both weak. The threshold reflectivity i?thr,si to 
overcome side loss is plotted in Fig. [TH for a pump with 
Iq = 10 15 W/cm 2 . A small reflectivity produces a large 
LW in the lower-right corner where damping is weak, thus 
allowing trapping to more easily occur. 

We now consider the specific plasma conditions men- 
tioned above as typical for SRBS on NIF ignition ex- 
periments, namely n e /n cr = 0.1 and T e = 2 keV 42 . 
The phase-matched SRBS modes have Ai = 553 nm and 
feAue = 0.297. The calculated backscatter gain rate is 
significant in both the CH ablator and He gas fill. The 
material affects a trapping assessment only via collisions. 
From Fig. QU for Z cS = 4 we find SN c /5N sl = 0.24, so 
we just consider side loss. The 3D side loss detrapping 
rate is z/a,si = 13.9/ps, or a time of l/^d.sl = 0.072 ps. 
The side loss threshold is SN s i — 2.7 x 10~ 3 , the critical 
intensity is I cr = 1.96 x 10 16 W/cm 2 , and the threshold 
reflectivity is fl t hr, s l = (5.28 x 10 13 W cm' 2 /I ) 2 . 

Finally, we show that smoothing by spectral dispersion 
(SSD) 33 is not likely to reduce trapping effects in SRBS 
on NIF. Recent experiments have utilized A/i = 45 
GHz of SSD bandwidth in the fundamental, 1054 nm 
laser light. After frequency-tripling, this corresponds to 
a speckle lifetime of t ss d = l/(3A/i) = 7.4 ps. For 
the reference SRBS conditions discussed above, t ss d ~ 
lOO/^si- Thus SSD is much less effective at detrapping 
than side loss. Moreover, a Langmuir wave overcomes 
SSD detrapping (tb < t ss d) for a very low amplitude 
of SN = 2.5 x 10~ 7 , or a reflectivity of (4.9 x 10 9 W 
cm- 2 // ) 2 . 



FIG. 12. (Color online.) Critical amplitude 8N cr from Eq. 
(|34|l for Z c ff=4. The side loss detrapping rate is faster than 
the collisional one for SN > SN cr . Blue curves are for SN cr = 
3 x 1(T 5 , 1 x 1(T 4 , 3 x 1CT 4 , 1 x 1(T 3 . 



The "critical intensity" I cr is introduced for convenience. 
Since trapping effects become significant for SN s» SN{ 
where SNi is the detrapping threshold for process i, we 
define the threshold reflectivity i?thr i for which SN = 
SN,: 



-Rthr.i = 



^fSNt 
Jo 



(40) 



To illustrate the threshold reflectivity, we consider 
SRBS of the phase-matched natural LW. The critical in- 
tensity I cr is plotted in Fig. Q21 We use x from Eq. (|T7|) 



VII. CONCLUSIONS AND FUTURE PROSPECTS 

This paper presented a framework for estimating when 
electron trapping nonlincarity becomes important in 
Langmuir-wave dynamics. Detrapping by convective 
loss in the longitudinal and transverse directions were 
discussed, as well as detrapping by Coulomb collisions 
(electron-electron and electron- ion) . 2D-2V simulations 
with the Vlasov code loki quantified trapping effects in 
driven LWs with finite transverse profiles, and showed 
they increase with the side loss bounce number as the 
transverse width increases. These runs showed trapping 
has little importance for kXoe = 0.7, which is above the 
loss of resonance value of 0.53. We explained this in terms 
of Eq. (flB)) , and related it to the lack of kinetic inflation 
in non-resonant Raman scattering. 

We applied our results to LWs driven by SRBS in NIF- 
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FIG. 13. Critical intensity I cr from Eq. (J39} in W/cm 2 re- 
lating laser intensity Io and reflectivity R to LW amplitude: 
6N = (Io/I cr )R 1/2 . 
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FIG. 14. Threshold reflectivity R tia , s l oc from Eq. ((40 
to overcome side loss for Io = 10 15 W/cm 2 . 



relevant conditions. An F — 8 intense laser speckle was 
taken as the LW size in order to estimate side loss. Side 
loss is generally more effective at detrapping electrons 
than collisions in speckles, although this is not the case 
for wider LWs. Linear response at the locally resonant 
SRBS scattered wavelength allowed us to obtain an lo- 
cal reflectivity needed for trapping to overcome side loss. 
This gives rather small values for plasma conditions from 
which SRBS is thought to originate in current NIF ex- 



periments. Moreover, it is the speckle, and not the lower 
beam average, intensity that matters. Preliminary as- 
sessment of pF3D simulations^ of NIF targets indicates 
a significant fraction of SRBS-generated LWs are above 
our trapping threshold. Future work will assess this, and 
attempt to incorporate trapping effects into enveloped 
propagation codes like pF3D. 
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Appendix A: Derivation of collisional results 



We restate our collision operator from Eq. (|23l) : 

-% = (i + z cS )u- 3 d,[(i-n 2 W] 

Vq ot 

+2u- 2 d u (f + u- 1 d u f). (Al) 

The second term describes collisions of tail electrons off 
bulk electrons, and is valid for u p ^> 1. The parallel (u x ) 
and perpendicular (u±) velocities are given by u x = fiu 
andw± = [l-/i 2 ] 1 / 2 uwithu G [0, oo] and fi e [-1,1]- We 
compute iVt r , the number of electrons initially trapped, 
that remain so up to time t. That is, once an electron 
leaves the trapping region its coherent bounce motion 
stops, even if it re-enters the trapping region later. The 
trapping region extends from u x = u p ± Ut r /2 and over 
all u±. 

Changing variables from to {u^,u x ) gives 



v dt ~ [Dxxd ^ 



+ D x d Ux 
D±d u± 



D ±± d\]f. (A2) 



The D's are straightforward to work out, and we do not 
give them. We assume u tr is small, and order derivatives 
as d/du x oj l/u tr and d/du± ~ 1 (v± ~ vt £ in physical 
units). For sufficiently small u tr , the dominant term is 
C xx d 2 f/du 2 x . For u p > 1, this is valid if u tr <§; F(Z)/u p 
where F(Z) is a function of Z. With this approximation, 
the collision operator yields a ID diffusion equation: 



v at 



' dui ' 



)u\{\ + Z cS ) + 2u 



{u\ + U 2 ) 5 /2 



(A3) 



(A4) 



We solve this equation subject to the outflow boundary 
conditions f(u x = u±,uj_,t) — with u± = u p ± u tr /2 
the boundaries of the trapping region. The initial con- 
dition for the trapped distribution is / — /oexp[— (u\ + 
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Up)/2] (a Maxwellian with u x evaluated at u p ) inside the 
trapping region, and / = otherwise. The number of 
trapped electrons is 



iVtr = 27T 



du i 



Vtr = ^Vl / UU± U 



du x f. 



(A5) 



We choose f = e u l /2 / (2i:u tT ) so N tr (t = 0) = 1. / has 
the solution 



/= ^2 fn(u±,t)si 



sm mrw, 



fn 



ra=l,3,... 

2 

ir 2 nu tr 



cxp[-u 2 ± /2 - n 2 Dt]. 



(A6) 



(A7) 



w = (u x — U-)/ut r , t is given by Eq. (1261) . and D = 
D xx u~ 3 . The sum is over odd positive integers since the 
even terms vanish. The trapped fraction becomes 



8 f°° 

— Vn~ 2 / dxexp\-x- n 2 Dt\. (A8) 
*■ „ Jo 



a; = /2 is a dummy integration variable. The decay 
rate of mode n goes like n 2 , as is typical of diffusion 
problems. After a short time, the n = 1 term dominates. 
Retaining just this term, and evaluating 8/-7T 2 = 0.81, we 
find 



A tr » 0.81/, 



/ = 



da; i 



-w 



W = x + Dt, 

(1 + Z)2x(l + 2ax) + 2 



D 



(1 + 2axf/ 2 



(A9) 
(A10) 
(All) 
(A12) 



To alleviate notation, we replaced t with t, Z c g with Z, 
and define a = w~ 2 . 

The upshot is an implicit integral equation for t: 



I(t,a) 



0.81 



(A13) 



Numerically finding t reveals it is linear in a for a < 1. 
We thus write t = to + aii and expand for a -C 1: 

J « J(t, 0) + ahd t I(t, 0) + a<V(t, 0) + 0(a 2 ) = 6 (A14) 

We choose to such that I (to, 0) = b. We find 

l(t,0) = f dxexp[-W ], (A15) 
VK = x + 2t [(1 + Z)x + 1] . (A16) 
Performing the integral gives an implicit equation for t: 

„-2t 

(A17) 



1 + 2(1 + Z)t 

This gives an exact formula for N tr in the limit u p — > oo : 
0.81e- 2 * 



JVtr(*) 



I + 2(1 + Z)t 



u v — >• oo. 



(A18) 



The temporal decay of A tr is thus not strictly exponen- 
tial. This formula reflects the different mathematical 
character of parallel dynamics from electron-electron col- 



lisions 



and pitch-angle scattering from collisions 



with all species ~ (1 + Z)t. Eq. (|A17I) is transcendental, 
and can be "solved" in terms of the Lambert W function. 
We are interested in cases where t < 1, so we Taylor ex- 
pand e~ 2t to order t 2 and obtain 



tn 



1 



Y +[Y 2 - 2 + 2b] 



1/2 ' 



Y = l + b(l + Z). (A19) 



This formula is valid (to real) for b above bo(Z). For 
Z = we have &o = 5 1 / 2 — 2 ps 0.236, and &o decreases 
with Z. For our choice of N tr = 1/2, b = 0.617 > 6 
for all Z. We have used the quadratic formula in a form 
that demonstrates the large- Z limit more clearly, which 
to leading order in Y is 



to 



1-6 

2F ' 



Y > 1. 



With iVtr = 1/2 this becomes 

0.31 



to 



2.62 



(A20) 



(A21) 



This form is accurate to within 10% for all Z > 0. The 
correction for finite a is 

h= -[d a I/d t I]\ t=kua=0 (A22) 

/°°exp[-Wo]a n Go 

- -to— Too ; ... i „ — ■ (A23) 

J Q exp[-W ]G 



The result is 



ti — to 



ll + 6Z+10i (l + Z) 



(1 + 2i (l + Z))(2 + Z + 2t (l + z)) ' 
Using Nt T = 1/2 and our approximate form for to , 
1.15 + 5.70Z- 1 + 6.09Z- 2 



Z + 7.23 + 16.3Z- 1 + U.7Z- 2 ' 



(A24) 



(A25) 
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